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INTRODUCTION 

In this paper we study the behavior of a queueing 
system which arises in the study of certain communication 
networks. Specifically we study a queueing phenomenon which 
arises with the SENET network, as described by Coviello 
and Vena (1975) or Barbacci and Oakley (1976) . This network 
allows for both voice and data messages to be transmitted 
over the same channels by using a special type of integrated 
circuit and packet-switched multiplexor structure. The two 
classes of traffic have substantially different performance 
requirements. Voice messages tend to possess great redundancy, 
and hence not to be sensitive to channel error rates, while 
data is very sensitive to channel error, having essentially 
no redundancy. Voice messages on the other hand have critical 
timing requirements and cannot be queued, while data is 
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relatively insensitive to timing and can be queued. Addition- 
ally, voice messages tend to be very long relative to data 
messages which can be broken up into small packets. These 
special requirements have led to the following queueing net- 
work. A node of the network consists of c + v channels or 
servers. The voice messages are assigned to v channels and 
do not queue. Thus the voice messages operate as a loss 
system. Data messages may use c channels exclusively and 
any unused voice channels; however, voice preempts data using 
voice channels. Data messages are queued if necessary. 

Typical performance measures that one may wish to calculate 
include the loss rate of voice traffic and the mean data 
queue length. 

We make standard probabilistic assumptions. Specifically, 
we assume voice traffic arrives according to a Poisson (A) 
process and each voice message has an independent exponential ( y) 
service time. Data messages are assumed to have independent 
exponential (n ) service times and arrive according to a Poisson(6) 
process. With these assumptions voice is an M/M/v/v loss 
system, and data is an M/M/S system where S=c+v-V(t) 
with V(t) = number of voice messages in service. The 
stochastic process { (X (t) , V( t) ) , t 0} is Markov with 
state space x {0,1,..., v} where X(t) = data system size 

at time t. One can easily write the Kolmogorov forward 
equations appropriate for this system; however, these equations 
do not yield a closed form solution. To describe this system 
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one must either numerically solve the forward equations or 
introduce approximations . 

This system has been studied previously by a number 
of researchers including Halfin and Segal (1972), Halfin (1972), 
Fischer and Harris (1976), Bhat and Fischer (1976), Fischer 
(1977) , Chang (1977) , and Gaver and Lehoczky (1979a, b) . The 
last two papers introduce a "fluid flow" and a diffusion 
approximation and derive explicit formulas for data queue 
behavior. These papers focus on the important case in which 
= 6/n > c. In such a situation the data messages must 
have access to voice channels for the system to be stable. 

4 

Furthermore, it was assumed that n/y was large, say 10 . 

Under these circumstances the data flow could be treated 
deterministically. Suppose we define = X/y and 

q = (p^/vI)/^j_Q p^/j I , the Erlang B blocking probability. 

The total traffic intensity on the c + v channels is given 
by P^ + p^(l-q) , or we could define p = (p^ + p^(l-q) ) / (c+v) . 
A heavy traffic approximation can be derived for this case 
p 1 . Such an approximation was derived in Gaver and 
Lehoczky (1979b) assuming n/p was large; a Wiener process 
with reflecting boundary was fo\and appropriate. In this 
paper we derive a heavy traffic approximation for the system 
without the fluid flow assumption that n/y is large. 

The methodology is drawn heavily from the approach of Burman 
(1979) . In this approach one characterizes a Markov process 



3 



by its infinitesimal generator. One next suitably normalizes 
the process so that the generator converges to a limiting 
infinitesimal generator (in this case to that of a reflected 
Brownian motion) . This convergence allows the conclusion that 
the finite dimensional distributions of the normalized Markov 
process converge. The diffusion approximation consists of 
treating the actual process through its limiting behavior. 

The details are somewhat complicated by the presence of a 
boundary . 
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2 . 



Let { (X(t) ,V(t) ) , t 0} be a bivariate Markov 
process with state space S = x {o , 1 , . . . , . Here 
{V(t), t 0} is marginally an M/M/v/v loss system with 
arrival rate X and service rate y- Conditional on V(t), 
{X(t), t 0} is an M/M/(c + v - V(t)) queueing system with 
arrival rate 6 and service rate n . We say that the V 
process subordinates the X process. We let 




( 2 . 1 ) 



the infinitesimal generator of the V process. 

The generator of the (X,V) process is given by 



Af (x, k) 



Qf(x,k) + (5{f(x+l,k) - f(x,k)) 

+ n(c+v-k) (f (x-l,k) -f (x,k) ) 
if X > c+v-k 



( 2 . 2 ) 



Of(x,k) +6(f(x+l,k) - f(x,k)) 

+ nx(f(x-l,k) - f(x,k)) 

if X < c+v-k 
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for f:S -> R continuous where 



Qf(x,k) = p^f(x,k+l) - (k+p^) f(x,k) + kf(x,k-l) 

V > k > 0 



(2.3) 



and f(x,-l) = f(x,v+l) = 0. Clearly Qf(x) = 0, that is Q 
annihilates functions of x alone. We next normalize the 
(X,V) process by defining ^^(t) = X(nt)//n and = V(nt). 

One can calculate the generator of the Markov process 
{(X^(t), t 0} having state space 

= { 0 , l//'n , 2//n , . . . } x { 0 , 1 , . . . , v) to be 



nQf (x,k) +6n(f(x+ l//li,k) - f (x,k) ) 



+ n^(f (x - l//“ n,k) 



if X > 



c + V - 



f (x,k) ) 
k 






A^f (x,k) 



nQf(x,k) +6n(f(x + l//n,k) - f(x,k)) 



(2.4) 



+ nn / nx (f(x + l//n,k) - f(x,k)) 



if X = 0, l//n,..., (c+v-k) //li’ . 



We assume f(x,k) has three bounded derivatives in x 
for each fixed k. With this assumption one can expand terms 
in (2.4) in a Taylor series and rewrite as 
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r 



nQf(x,k) + n^/^f^(x,k) (6-n (c+v-k) ) 



+ 2 f (x,k) (6 +n (c+y-k) ) 



Aj^f(x,k) 



+ 0(n-l/2) 



if X > (c+v-k) //n 



nQf(x,k) + (x, k) (6 -ri/ nx) 






(2.5) 



if X = 0, l//n, ... , (c+v-k) //n 



with f^(x,k) = — f(x,k) and f^^(x,k) = — ^ f (x,k) . 

3 X 

We ultimately wish to prove that the finite dimensional 
distributions of {X^(t), t 0} converge to those of a 
Wiener process with reflecting barrier at the origin. This 
can be restated in terms of semi-groups. We let t >_ 0} 

be the semi-group of operators associated with 
{(X^(t), Vj^(t))r t 0} and t >_ 0} be that associated 

with a Wiener process having reflecting barrier at 0. Let 
g be a continuous function g:R' ->■ R' . Knowledge of the 
semi-group is equivalent to knowledge of the transition 
functions by taking a sequence of g's which approximate 
indicator functions. We wish to prove I g(x,k) - g(x)| ^0 

as n ^ for all (x,k) . Here g(x,k) = E(g(X^) |Xj^(0) = x, 

Vj^(O) = k) . The presence of the variable k prevents this 



7 



from being done directly. The method we use is to construct 

oo 

a convenient sequence of functions ^^n^n=l '^hich converge 
in some sense to g. We write 



iT^g 



- T 



gii < 

:^"n — 



'Vn 



- VI 



n 



iT^g - T^g li 
I t^n"n 



iT^g - (T'”g) li 
' t^n ^ t^^n"n 



( 2 . 6 ) 



where || || refers to the sup norm over . Both and 

n n t 

oo 

are contraction semigroups. 

oo oo , 

<(T^g)^>^_j^ is the sequence of functions constructed 

oo 

from T^g* goal is to show that each of the three terms 

on the right side of (2.6) converges to 0. The first and 

second terms can be handled similarly. For any function g, 

oo 

we must guarantee that the constructed ^?n^n=l sequence 

II M ,, oo oo ,, 

satisfies g - g -^0. It will follow that (T.g) - T.g 0. 

^n n t nt n 

Moreover, since t >_ 0} is a contraction semi-group 

liT^g - £ ilg “ ^n^^n '^hich also converges to 0. The 

oo 

sequence ^^n^n-1 chosen in such a way that the third 

term converges to 0 . 

We focus on a convergence determining class of functions 
g, those which are bounded and have three bounded derivatives. 

For such a function g(x) we define 



gj^(x,k) = g(x) + 



-i— u(x,k) + ^ v(x,k) 

/IFT 



(2.7) 
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where u and v have two bounded derivatives in x for 
each fixed k. The functions u and v will be determined 
explicitly later and are chosen to control the third term in 
(2.6). Clearly when is defined by (2.7), = 0(n~^^^) 

and therefore converges to 0 as required. 

One can apply the generator to g^ to derive 



nOq(x) + n^-^^[Qu(x,k) + g ' (x) (6 -n (c+v-k) ) ] 
+ [Qv(x,k) + u (x,k) (6 -n (c+v-k) ) 

+ I g" (X) (6 -n (c+v-k) ) ] + 0(n"^/^) 
if X > ( c+v-k) //^n 



Angn(x.k) 






nOg (x) + (Qu{x,k) + g'(x)(5- nn^^^x) ] 

+ [Qv{x,k) + Uj^ (x,k) {6 -nn^^^x) + j g” (x) (6 +nn^^^x) ] 



+ 0(n-l/2) 



if X = 0, l//n,..., (c+v-k) //ri 



g 

where u^(x,k) = u(x,k) . Recall that Q annihilates functions 

of X alone, thus nQg(x) =0. We want to have A^g^^(x,k) 

converge to a finite limit and to have that limit be inde- 

1/2 

pendent of k. For this to occur, the n ' term must be 
controlled and the functions u and v must be chosen in 
such a way as to eliminate the variable k. 
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1/2 

The n coefficient in (2.8) can be rewritten by 

adding and subtracting 

V 

I TT, g' (x) (6 -n (c+v-k) ) = -r) (c+v) (1-p) g' (x) . 

k=0 ^ 

We next pick u(x,k) to be a solution of 

Ou(x,k) = -(g'(x) n (p^- (c+v-k) ) + g' (x) n (c+v) (1-p)) 
= -g' (x) n(k-p^d-q) ) (2.9) 

When u(x,k) is any solution of (2.9), the coefficient of 
1/2 

the n term in (2.8) becomes 



-g' (x) n (c+v) (1-p) 



if 



X > HlZli 

/ n 



g' (x) n((c+v)p - n^^^x - k) 



if 



0 < X < 



c+v-k 
/ n 



Equation (2.9) can be solved explicitly. Define 
a, = -g' (x) ri(k-P (l-q))/M, so that (2.9) can be written as 

K V 

-p^(u(x, k) -u (x,k-l ) ) - (k-l) (u (x,k-l ) -u (x,k-2) ) = / k = l,...,v 

-v(u (x,v) -u(x, v-1) ) = a^ (2.10) 
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0 , 



Equation (2.10) has a solution since I^_q ~ 



<V>Lo 



where <tt, is the stationary distribution associated 

k 



with Q, or 7T, = (p /k'. ) / (T p^/i'.). The solution is 

K. V 1*“U V 

given by 



u(x,k) - u(x,k-l) = 



k-1 

y TT . a . 
i=0 ^ ^ 

Pv\-1 



-g' (x)nTj^_^ 



k-1 



where 



k 

T = I TT (1-p (1-q)) and T = 0. 

K i^O 1 V V 



Clearly 



u(x,k) = u(x,0) 



g' (x)n 
yp 

V 



k 

I 

i=l 



!izi 

^i-i 



1 < k < V 



( 2 . 11 ) 



where u(x,0) is arbitrary. We let u(x,0) 



I 9'(x) 



SO 



u(x,k) 



g' (X) 



yp 



k 

- I 

V i=l 



Vi 

^-1 



0 < k < V (2.12) 



For the choice of u specified by (2.12) we next wish 
to insure that the limiting generator is independent of the 
variable k. The function v is chosen to eliminate the 
dependence on k. The 0(1) term of (2.8) is given, for 
X > (c+v-k)//n, by 
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Qv(x,k) + g" (x) 



1 

2 



k T 

I (<S -n (c+v-k)) 



y p . ^ , 7T . , 

1=1 1-1 



+ J g" (x) (6 +n (c+v-k) ) 



= Qv(x,k) + H (x,k) . 



Let H(x) = TTj^H(x,k) and consider Qv(x,k) + (H (x,k) -H (x) ^+H (x) 

We now let v(x,k) be any solution of 



Qv(x,k) = -(H(x,k) - H(x)) . 



(2.13) 



Equation (2.13) has a one-parameter family of solutions, 
since “ H(x)) = 0. When v(x,k) is chosen 

to be any solution of (2.13), the 0(1) term of (2.3), 

for X >; (c+v-k) //n, will become H(x) and will therefore 
be independent of k. It remains to calculate H(x) . 



H (x) = g" (x) 



V 



I 



k=0 



yp 



k 

- I 

V i=l 



i-1 

^i-1 



(5 -n (c+v-k) ) 



+ -2 (6 +n (c+v-k)) 



= g"(x) 



6 - I 7T, (6 -n (c+v-k) ) I 



V 



k T 



yp 



V k=0 



i-1 



1=1 1-1 



(2.14) 



= g" (X) 



2 V 



6 - 



yp , ^ n k 
^■^v k=0 



I 7t, (k-p (1-q)) I 



k T._^ 



i=l ^i-1 



n (c+v) (i-p) y ^ y ^i-l 
yp 1 A "k . TT . , 

^^v k=0 1=1 1-1 
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The second term can be rewritten by interchanging the order 
of summation. The third term is O(l-p) . We find 



H(x) = g"(x) 



2 v-1 T. V 



6 - 






ypv i=0 ^i k=i+l 



= g" (x) 



2 v-1 T . 



6 - 



yp 



y — (T -T.) + O(l-p) 

^ ri TT . V 1 



V i=0 i 



with T =0 or 

V 



H(x) = g" (x) t] 



v-1 T7 

p, + I -i + O(l-p) 

d up . IT . ' K / 

^^v 1=0 1 



(2.15) 



For the functions u and v specified by (2.12) 
and (2.13), equation (2.8) can be rewritten as 



A q (x,k) 
n^n ' 



(1-p ) (c+v) ng' (x) 



T v-1 T7 

pj + — y — — + o(i-p) 

yp^ i£o ^i 



g"(x) + 0(n ^/^) 



for X > (cfv-k)//n 



= \ 



(2.16) 



n[(c+v)p -n^'^^x-k] g' (x) 



+ ng"(x) 



+ 



d yp 



v-1 T7 

I — + O(l-p) 

^ . TT . 



V i=0 ^i 



I/O k T. 

- (c+v-k-n '^^x) y ^ 

^Pv i=0 ^i 



t 0(n-l/2) 



for X < (c+v-k)//n 
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We now introduce the "heavy traffic approximation." 

In order for the generator to converge to a limiting 

- 1/2 

generator we must have 1-p = 0(n ' ) . Specifically, we 
assume p = = 1 - (0//n) for some 0 0 . In this case, 

n^/^(l-p) = e, and (2.16) becomes 






-0n ( c+v) g* (x) + n 



P . + 

d pp 



v-1 T?“* 

I — 

V i=0 ^i 



g"(x) + 0(n"^/^) 



A g (x,k) = 
n^n 



for X (c+v-k)//n 

n [ (c+v) p-n^'^^x-kjn^^^g* (x) 

+ ng"(x) 

- (c+v-k-n^'^^x) 





D + *^ 


v-1 

y 


t2] 
1 . 




yp 

L V 


L 

i=0 


TT . 

1 J 



k T. 

n y _i 

y P . j ^ i 



V 1=0 1 

for X < (c'+v-k)//n 



(2.17) 






We now define a limiting generator with domain 

consisting of all functions g having three bounded derivatives 
and g' ( 0) = 0 . Let 



A g(x) = -0n(c+v) g' (x) + n 



p j + — 

d yp 



v-1 T. 

I — 

V i=0 ^i 



2 -1 



g" (X) , 



X > 0 



(2.18) 
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is the generator of a Markov process which 
corresponds to a Wiener process with drift -0n(c+v), scale 



2n 




+ 




v-l 



I 

i=0 




/ 



- 1/2 

and a reflecting barrier at 0. The 0(n ) terms involve 

the first three derivatives of g which are bounded. It is 
clear from a direct comparison of (2.16) and (2.18) that 
|a g (x,k) - Agl^O as n->« for all x > 0 and k 
arbitrary. In addition, g' (0) =0 is necessary for the 
generator to converge at x = 0. Unfortunately even assuming 
g' (0) = 0, 

2 k T^ 

|Aj^g(0,k) - A^g(0)| ^ (c+v-k) ^ g"(0) as n ^ « 



rather than to 0. One needs a special argument to handle 
this lack of convergence at the boundary. 

We set out to prove the third term in (2.6) converges 
to 0. Standard semi-group results (see Burman, 1979, p. 33) 
give 



oo n 

(T.g) - T.g 

t^ n t^n 



= / 

0 






-Aw ) dS 
n n 



(2.19) 



where w = w(t,x) = T~g(x) . Recall that w^ = w +(l//n) u+(l/n) v 
with u and v defined by (2.12) and (2.13) with g replaced 
by w. It follows that 
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It^q ” (T q) II 

' t^n ^ t^^n"n 



t 

= 11/ Tj,.g((A_^w)^-A^W+A„w-A^W^)dS||^ 

< / l|T’J_g((A^w)_^-A^w)l|„dS + 11/ T^_g(A„w-A^Wj,)dS||j^ 

< / ll(A„w)^-A„w||^dS + 11/ T^.3(A„w-A„w^)dS||_^ . 



- 1/2 

The first term is clearly 0(n ' ) . It remains to 

- 1/2 

show that the second is 0(n ) as well. We have shown 

|Aw-Aw|=0(n ' ) except at the boundary where it 
oo n n 

is 0(1). We split the integral into two parts, for one 

of which the process is away from the boundary, and for the 

other, near the boundary. The integral away from the boundary 

- 1/2 

has an integrand which is 0(n ' ). The integral near the 

- 1/2 

boundary is also 0(n ) since under a heavy traffic 

assumption the process is rarely near the boundary. The 

details are merely summarized here; they are based on the 

ideas of Burman (1979) . 

Let I be the indicator function of 
on 

and I. be the indicator of 
In 
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We have 









<11/ (A„w - A„w„)I^„dsll„ + llA^w - A^w^ll„ 11/ T^VonlSlIn • 

“ 1 /2 — 1 /2 
The first term is 0(n ' ), since |a w - A w | = 0(n ' ) off 

n n 

the boundary. The factor IIa^w - A^w^H = 0(1), thus it remains 
to show that 



0 






This gives the total time in [0,t] spent near the boundary. 
We bound 



by first introducing a function h(x) not in the domain of 

A^. We let h(x) have bounded support, be infinitely 

differentiable and be given by h(x) = x for x near 0. 

One can construct h (x) using (2.7) and apply A to 

h to find 
n 



0 ( 1 ) 



if X > 



c+v-k 

/n 



Ah = 
n n 



( 2 . 20 ) 



( (c+v) p-n^^^x-k) + 0(1) if x < ^ 

/ n 



One has 



A - h = / T^A h dS 
t n n ^ S n n 



t t 

= / T„A h I dS + / T^A^h I dS. 
Q S n n J.n q S n n 



It follows that 



t"a h I. dS|| < ||T“h - h II + || / T“a h I, dS|| 

S n n On "n — " t n n"n "i S n n In "n 



„n. 



„n. 



< 2||h II + 0(1) 
— " n"n 



We have shown || / T^A h dA|| to be bounded in n. An 

S n n On "n 



application of (2.2) shows 



TgA^hj^Ip^dAll^ = n^/^n I (C+v) p - n^/^x - k + 0(1) I ||/ Tgl^^dsll^^ 



is bounded in n. It follows that ||/ "^s^on^^^^n" ^'^^) • 
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This finally concludes the argument which shows 

■ ^“^t^^nlln " 0(n“^/^), hence by (2.6) ||T^g - T~g||^ = 0(n"^/^) 

We have thus shown that the finite-dimensional distri- 
butions of the (X^(t) , process converge to those 

of a Wiener process with drift -9n(c+v) scale 



n 



( v-l T? 

p + _n_ I ^ 
V iio "i 



f 



and reflection at 0. The diffusion approximation treats 
X^(t) as though it were such a Wiener process. For instance, 
the limiting Wiener process has a stationary exponential dis- 
tribution with parameter 



0 (c+v) 
— 

I 

V i=0 



+ 



_n 

UP 



(TJ/tt.) 



This is a distribution for X(nt)//n and suggests X(t) 
will have a steady state distribution given approximately 
by an exponential with parameter 

(c+v) U-P)/(p<j + ^ 

The steady state mean data queue length would then be 



v-l T^ \ 
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E (X(t) ) 



( 2 . 21 ) 



Pd ^ 



2 

„ v-1 

n y _i 

yp . TT • 
1=0 1 



(c + v) (1-p) 



It is interesting to consider the special case c = 0, 

V = 1 where the two types of traffic use the same channel. 
Under heavy traffic p = p^ +p^/(l+p^), so P^j ~ 

The mean data queue length derived from the diffusion approxi- 
mation (2.21) will be 




The latter is the exact expression derived by Fisher (1978) 
for this case. The expression (2.21) represents a generaliza- 
tion of the results of Gaver and Lehoczky (1979b) . In this 
paper, a diffusion approximation is given based on the fluid 
flow assumption for the data. For this case the result is the 
same except that the scale is given by 



rather than 




v-1 

I 



i=0 




7T . 

1 



n 




+ 




v-1 



I 

i=0 




The results derived in this paper therefore definitely generalize 
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the results of Gaver and Lehoczky (1979b) since the variability 
in the data queue is now included. When n/y is large, the 
second term dominates, and the fluid flow approximation is 
satisfactory . 



provides a method for studying the dynamics of that process. 

For instance, suppose the data queue were at level x at 
time t where x is large. One might wish to study the time 
that elapses until the queue becomes empty. This is essentially 
the duration of the busy period under heavy traffic and 
corresponds to a first-passage time for a Wiener process. Let 
us denote it by T . Straightforward martingale arguments 
provide for its transform 



The Wiener process approximation for the X(t) process 



-sT 

E (e ^ 




( 2 . 22 ) 



where 



m = 0(c+v)n ~ n^'^^ (1-p) (c+v)n 




It is also easy to find the mean first-passage time 



E (T ) = x/m 

X 



(2.23) 
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One might also be interested in the area beneath the 



sample path until emptiness occurs, since this area represents 
the total time waited by all data customers involved in the 
busy period. If A represents this area, simple backward 
equation arguments give 



E(A^) 




where 



m and 



2 

a 



are given in (2.22). 



( 2 . 24 ) 
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